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^ , Abstract. It is a major unsolved problem as to whether unknot recogni- 

j— _ i tion — that is, testing whether a given closed loop in R 3 can be untangled to 
| ■ form a plain circle — has a polynomial time solution. In practice, trivial knots 

(which can be untangled) are typically easy to identify using fast simplification 
techniques, whereas non-trivial knots (which cannot be untangled) are more 
resistant to being conclusively identified as such. Here we present the first 
1 unknot recognition algorithm which is conclusive and, although exponential 

Otime in theory, exhibits a clear polynomial time behaviour in practice even 
for non-trivial knots. Specifically, through exhaustive testing over all 2977 
pH , non-trivial prime knots with < 12 crossings, the algorithm runs in every case 

by solving only a linear number of linear programs. The algorithm draws on 
techniques from normal surface theory and combinatorial optimisation, and 
we also adapt it to the important topological problems of 3-spherc recognition 
and prime decomposition. 
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■ 1. Introduction 

One of the most well-known computational problems in knot theory is unknot 
recognition: given a knot K in R 3 (i.e., a closed loop with no self- intersections), 
can it be deformed topologically (without passing through itself) into a trivial 
unknotted circle? If the answer is "yes" then K is called a trivial knot, or the 



X 



unknot (as in Figure 1(a) ); if the answer is "no" then K is called a non-trivial knot 
(as in Figure 1(b) I. This simple yes/no decision problem is deceptively complex: 
the best known algorithms require worst-case exponential time, and it is currently 
a major open problem as to whether a polynomial time solution is possible. 

The input to this problem is typically a knot diagram, i.e., a piecewise-linear 
projection of the knot K onto the plane in which line segments "cross" over or 
under one another, as seen in Figure [TJ The input size is typically measured by the 
number of crossings c (note that a c-crossing knot can be deformed into a diagram 
with O(c) line segments). 

Only little is known about the computational complexity of this problem. Un- 
knot recognition is known to lie in NP [23], and also in co-NP if the generalised 
Riemann hypothesis holds [30]. Haken's original algorithm from the 1960s [21] has 
been improved upon by many authors, but the best derivatives still have worst-case 
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(a) Some trivial knots 



(b) Some non-trivial knots 



FIGURE 1 . Examples of knots in R 3 



exponential time. There are alternative solutions, such as Dynnikov's grid simplifi- 
cation algorithm [T5], but these are likewise exponential time or worsen Neverthe- 
less, if we assume the generalised Riemann hypothesis and that NP ^ co-NP then 
we have reason to believe that unknot recognition should not be NP-complete, and 
nowadays there is increasing discussion as to whether a polynomial time solution 
might indeed exist [HI [22] . 

For inputs that are trivial (i.e., topologically equivalent to the unknot), solving 
unknot recognition appears to be easy in practice. There are widely-available sim- 
plification tools that attempt to "reduce" the input to a smaller representation of 
the same topological knot in polynomial time [3J [TT] , and if these tools can reduce 
the input all the way down to a circle with no crossings then the problem is solved. 
Experimentation suggests that, given modern simplification tools, it is extremely 
difficult to find "pathological" representations of the unknot that do not simplify 
in this way [3] [10] . 

For input knots that are non-trivial (i.e., inequivalent to the unknot), the situa- 
tion is more difficult. Here our simplification tools cannot help: they might reduce 
the input somewhat, but we still need to prove that the resulting knot cannot be 
completely untangled. There are many computable knot invariants that can assist 
with this task |2j [32] , but all known invariants either come with exponential time 
algorithms or might lead to inconclusive results (or both). 

In this sense, obtaining a "no" answer — that is, proving a knot to be non-trivial — 
is the more difficult task for unknot recognition. Indeed, this mirrors the earlier 
complexity results, in which unknot recognition has been known to lie in NP for 
over a decade, whereas the corresponding co-NP result is much newer and relies 
on an unproven mathematical conjecture. 

Our contribution is a new algorithm for unknot recognition, drawing on tech- 
niques from topology, combinatorial geometry and linear optimisation, that exper- 
imentally behaves like a polynomial time algorithm even when the input knot is 
non-trivial. In particular, when run over the Knotlnfo database of all 2977 prime 
knots with < 12 crossings [14], it proves each of them to be non-trivial by solving 
a linear number of linear programs. Combined with the aforementioned simplifica- 
tion tools, this gives the first algorithm for unknot recognition that guarantees a 
conclusive result, and that in practice exhibits polynomial-time behaviour. 

More precisely: we begin with a c-crossing diagram of the input knot K, and 
construct an intermediate 3-dimensional triangulation T with n tetrahedra, where 
n G 0(c). The algorithm then recursively builds and solves a sequence of linear 



Musick posted a preprint in 2011 that claims a polynomial time solution, and made a sub- 
stantial rewrite in 2012 |35| . The status of his claim remains unclear at present. 
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programs (each in polynomial time), essentially searching for a certificate of unknot- 
tedness. The key observation is that for all 2977 input knots, the total number of 
linear programs solved is approximately 8n. Combining these factors, we find that 
for all 2977 input knots the total running time grows at a rate that is polynomial 
in the input size. 

We emphasise that, like all known algorithms for unknot recognition to date, the 
theoretical running time remains worst-case exponential. For our new algorithm, 
the exponential running time stems from the total number of linear programs, 
which in theory could grow as large as 0(3 n • n 2 ). Nevertheless, the experimental 
performance is far better: the input data set of 2977 knots is an exhaustive list (not 
cherry-picked), yet each and every case exhibits this polynomial time behaviour 
(~ 8n linear programs in total). 

We also emphasise that this polynomial time behaviour is measured purely by 
experiment. We do not prove average-case, smoothed or generic complexity re- 
sults; although highly desirable, such results are sorely lacking in the realm of 
3-dimensional triangulations. There are good reasons for this scarcity of results, 
which we discuss further in Section |6l 

The new algorithm, described in Section [3] is based on normal surface theory, 
a powerful algorithmic toolkit at the heart of Haken's original unknot recognition 
algorithm [21) . Typical normal surface algorithms are based on combinatorial enu- 
meration, but we replace this with combinatorial optimisation following early ideas 
of Casson and Jaco et al. [26] . Essentially, we minimise the genus of a surface in R 3 
that the knot can bound: this genus is zero if and only if the input knot is trivial. 

Optimisation-based approaches have been attempted before, but none have ex- 
hibited polynomial time behaviour. The key difficulty is that we must optimise a 
linear objective function (which measures genus) over a non-convex domain (which 
encodes potential surfaces). In previous attempts: 

• Casson and Jaco et al. split the domain into an exponential number of 
convex pieces and run a linear program over each [26 . This gives a useful 
upper bound on the running time of 0(3™ xpoly(n)). However, the resulting 
algorithm is unsuitable because for non-trivial input knots the running time 
also has a lower bound of 0(3™ x poly(n)), whereas alternative algorithms 
have a slower growth rate in practice |10) . 

• The authors express this optimisation as a high-dimensional integer pro- 
gram 13J . The resulting integer programs yield effective bounds on the 
minimum genus, which are useful for related problems such as knot genus 
and crosscap number. However, these integer programs are extremely diffi- 
cult to solve exactly, which makes them unsuitable for unknot recognition. 

In our algorithm, we overcome this non-convexity issue using branch-and-bound 
techniques, where we run a sequence of linear programs over a combinatorial search 
tree. The broad layout of this search tree is inspired by recent normal surface 
enumeration algorithms [12], and its detailed structure is carefully constructed — 
and dynamically reorganised as the algorithm runs — to effectively discover and 
exploit additional constraints as they appear within the input triangulation. The 
triangulation itself is also modified to further streamline and constrain this search 
tree. Both this modification of the triangulation and the dynamic reorganisation 
of the search tree are crucial for the experimental polynomial-time behaviour that 
we observe in practice. 
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In addition to unknot recognition, we adapt our new algorithm to the related 
topological problems of 3-sphere recognition and prime decomposition of 3-manifolds. 
We discuss these applications in Section 01 and in particular we run the resulting 
new 3-sphere recognition algorithm over a comprehensive collection of 150 inputs 
(all of which are non-trivial homology spheres). Here the new algorithm shows 
an exponential time profile, but nevertheless runs orders of magnitude faster than 
prior state-of-the-art alternatives (9J [12] . 



2. Preliminaries 

Here we give a very brief overview of knots, triangulations and normal surfaces. 
For more details on the role of normal surfaces in unknot recognition and related 
topological problems, see the excellent summary by Hass et al. [23] , 

We consider a knot to be a piecewise linear simple closed curve in IR 3 , formed from 
a finite number of line segments. Two knots K, K' are considered equivalent if one 
can be continuously deformed into the other without introducing self-intersections. 
Any knot equivalent to a simple planar polygon is said to be trivial, or the unknot; 
all other knots are said to be non-trivial. Again, see Figured] for examples. 

A knot diagram is a projection of a knot into the plane with only finitely many 
multiple points, each of which is a double point at which two "strands" of the 
knot cross transversely, one "passing over" the other. These double points are 
called crossings: the four knot diagrams in Figure [T] have 0, 2, 3 and 4 crossings 
respectively. Alternatively, knot diagrams can be described as annotated 4-valent 
planar multigraphs; see [23 for the details. A knot diagram with c crossings can 
(up to knot equivalence) be described in 0(c) space (see [25] for some examples of 
encoding schemes), and in this paper we treat knot diagrams as the usual means 
by which knots are presented as input. 

By adding a point at infinity, we can extend the ambient space from R 3 to the 
topological 3-sphere S 3 = M. 3 U {oo}. For any knot K, we can then remove an open 
regular neighbourhood of K from S 3 (essentially "drilling out" the knot from S 3 ); 
this yields a 3-manifold with torus boundary called the knot complement K. Given 
a knot diagram with c crossings, Hass et al. show how to construct a triangulation 
of K with 0(c) tetrahedra in 0(c log c) time [23] Lemma 7.2]. 

Although the Hass et al. construction produces a simplicial complex, computa- 
tional topologists often work with generalised triangulations, which are more flexible 
and often significantly smaller. A generalised triangulation begins with n abstract 
tetrahedra, and affinely identifies (or "glues") some or all of their 4n triangular 
faces in pairs. Two different faces of the same tetrahedron may be glued together; 
moreover, as a consequence of the face gluings we may find that multiple edges of 
the same tetrahedron become identified together, and likewise with vertices. Of 
particular importance are one-vertex triangulations, in which all An tetrahedron 
vertices become identified as a single point. 

Essentially, these are the 3-dimensional analogues of well-known constructions 
in two dimensions, such as the two-triangle torus and the two-triangle Klein bottle 
shown in Figure [2] (the different arrowheads indicate how edges are glued together). 
Both of these are one- vertex triangulations. 

Unless otherwise specified, all triangulations in this paper are (3-dimensional) 
generalised triangulations. 
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(a) A torus (b) A Klein bottle 

Figure 2. Triangulated surfaces in two dimensions 



Figure 3. Normal triangles and quadrilaterals 

& # 4> # 

Figure 4. The seven normal disc types 

If a triangulation T represents some underlying 3- manifold M. , then the tetrahe- 
dron faces of T that are not glued to anything form a triangulated surface (possibly 
empty, possibly disconnected) which we call the boundary of T, denoted dT; this 
in turn represents the manifold boundary dM.. 

A normal surface in T is a surface 5* that is properly embedded (i.e., embedded 
so that dS C dT), and which meets each tetrahedron of T in a (possibly empty) 
collection of curvilinear triangles and quadrilaterals, as illustrated in Figure [3] In 
each tetrahedron these triangles and quadrilaterals are classified into seven types ac- 
cording to which edges of the tetrahedron they meet, as illustrated in Figure 2J four 
triangle types (each "truncating" one of the four vertices), and three quadrilateral 
types (each separating the four vertices into two pairs). 

Any normal surface S in an n-tetrahedron triangulation T can be described 
by a vector of 7n non-negative integers that counts the number of triangles and 
quadrilaterals of each type in each tetrahedron; we denote this by v(S') £ l7 n . The 
individual coordinates of v(S) that count triangles and quadrilaterals are referred to 
as triangle and quadrilateral coordinates respectively. This vector uniquely identifies 
the surface (up to a certain class of isotopy). More generally, an integer vector 
x e R 7n represents a normal surface if and only if: 

(1) x>0; 

(2) Ax. = 0, where A is a matrix of < 6n linear matching equations derived 
from the specific triangulation 7~; 

(3) x satisfies the quadrilateral constraints, which require that for each tetra- 
hedron at most one of the three corresponding quadrilateral coordinates is 
non-zero. 

In essence, the matching equations ensure that normal triangles and quadrilat- 
erals can be matched together across adjacent tetrahedra, and the quadrilateral 
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Figure 5. A disc bounded by a trivial knot 

constraints are required because any two quadrilaterals of different types in the 
same tetrahedron must intersect (which an embedded surface cannot have). Any 
vector in R 7n (real or integer) that satisfies all three of these conditions is called 
admissible. 

Haken's original unknot algorithm is based on the observation that any trivial 
knot must bound an embedded disc in R 3 (see Figure [5]). In the knot complement 
K, this corresponds to a properly embedded disc in K that meets the boundary 
dK in a non-trivial curve (i.e., a curve that does not bound a disc in the torus dK). 
Moreover, we have: 

Theorem 1 (Haken) . Let K be a knot and let T be a triangulation of the comple- 
ment K . Then K is trivial if and only if T contains a normal disc whose boundary 
is a non-trivial curve in dT. 

Let V be a vertex of some triangulation T ■ Then the link of V is the frontier of 
a small regular neighbourhood of V. If T triangulates a 3-manifold then this link 
is a disc (if V lies in the boundary dT) or a sphere (if V lies in the interior of T) ■ 
The link can be presented as a normal surface built from triangles only, and if T 
is a one- vertex triangulation then this link contains precisely one triangle of each 
type. 

We finish this section with the well-known observation that, in a fixed n-tetrahedron 
triangulation T, the Euler characteristic of an arbitrary normal surface S can be 
expressed as a homogeneous linear function \ : R 7 ™ — > R. There are many choices 
for such a function (since we arc only interested in its values on the admissible 
points of R 7 ™); see [13] for one such formulation. 

3. The algorithm 

Broadly speaking, the new algorithm operates as follows. 

Algorithm 2 (Unknot recognition). To test if an input knot K is trivial, given a 
knot diagram of K: 

(1) Build a triangulation of the knot complement K. 

(2) Make this into a one-vertex triangulation T using subcomplex expansion 
and crushing (described in Sections \3.1\ and \3.2\l , without increasing the 
number of tetrahedra. Let the number of tetrahedra in T be n. 

(3) Search for a connected normal surface S in T which is not the vertex link 
and has positive Euler characteristic. The search procedure (described in 
Section \3.S\) runs a series of linear feasibility tests over a combinatorial 
search tree that we build dynamically, and is motivated by more general 
branch- and-bound techniques. 

(4) If no such S exists, then the knot K is non-trivial. 

If S is a disc whose boundary is a non-trivial curve in dT, then the knot 
K is trivial. 
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FIGURE 6. Cutting along and then crushing a surface 




Figure 7. Pieces obtained during crushing 

Otherwise we modify T by crushing the surface S to obtain a new trian- 
gulation T' of K with fewer than n tetrahedra ( as detailed in Sections \3.1\ 
and \3.4\ l, and return to s£ep[H using this new triangulation T' instead. 

The remainder of this section is devoted to the details for the individual steps, 
which draw on a variety of techniques from combinatorial optimisation, computa- 
tional topology and normal surface theory. Moreover, we prove: 

Theorem 3 (Correctness). Algorithm^ is correct in determining whether K is 
trivial or non-trivial. 

Theorem 4 (Running time). Let c be the number of crossings in the input knot 
diagram for Algorithm^ With the exception of step [3] (where we search for the 
normal surface S ), every step of the algorithm runs in time polynomial in c; in 
contrast, step\3\ runs in time 0(3™ x poly(n)), where n is the number of tetrahedra 
in the triangulation T . Every step (including step\&$ is repeated at most O(c) times. 

As noted in the introduction, this makes the algorithm exponential time over- 
all. However, the search procedure in step [3] is sufficiently powerful to exhibit 
polynomial-time behaviour in practice (as quantified experimentally in Section [5]), 
and Theorem [4] then ensures that this translates to polynomial-time behaviour for 
the algorithm as a whole. 

We complete the proofs of Theorems and 0] in Section 13.51 (after providing the 
details of the algorithm). 

3.1. The Jaco-Rubinstein crushing process. Steps [2] and Q] make use of the 
Jaco- Rubinstein crushing process [27j , which modifies a triangulation by "destruc- 
tively" crushing a normal surface within it. Here we outline and analyse this pro- 
cess. 

The Jaco-Rubinstein technique works as follows. Given a normal surface S within 
a triangulation T, we cut T open along 5* and then crush all of the triangles and 
quadrilaterals of S (which now appear twice each on the boundary) to points. See 
Figure[6]for an illustration. This results in a new cell complex built from tetrahedra, 
footballs and/or pillows, as illustrated in Figure [7] We then flatten each football 
to an edge and each pillow to a triangular face, as shown in Figure [SJ resulting in 
a new triangulation T' . This new triangulation might be disconnected, and might 
represent a different 3-manifold from T (or indeed might not represent a 3-manifold 
at all). 

The tetrahedra in T' correspond precisely to those tetrahedra in T that did not 
contain any quadrilaterals of S. In particular, unless S is a collection of vertex 
links, the number of tetrahedra in T 7 will be strictly smaller than in T ■ 
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Figure 8. Flattening footballs and pillows 

A key result of Jaco and Rubinstein is that, if S is a sphere or a disc and the 
underlying 3-manifold is sufficiently simple, then any topological changes that result 
from this crushing are of limited types, and can be easily detected. In our setting, 
we use this to prove: 

Lemma 5. Given a triangulation T of a knot complement K with n > 1 tetrahedra 
and a connected normal surface S in T which is not a vertex link and which has 
positive Euler characteristic, we can crush S using the Jaco- Rubinstein technique 
and convert the result into a new triangulation T of K with strictly fewer than n 
tetrahedra. Moreover, this process can be carried out in 0(n) time. 

Proof. We briefly outline the key ideas here; see the appendix for the full proof. 

First we observe that S must be a sphere or a disc (as knot complements do not 
contain embedded projective planes). 

By results of Jaco and Rubinstein [27] , using the fact that T represents a knot 
complement, we show that any topological changes are limited to (i) cutting open 
the complement of a trivial knot to make a 3-ball; and (ii) adding or removing 
connected components that are 3-balls or 3-spheres. Therefore, if there is some 
component remaining with torus boundary then this becomes our new triangulation 
T 1 ; if not then we know that the input knot is trivial. 

Since S is connected but not a vertex link, some tetrahedron of T contains 
quadrilaterals of S, and so the crushing operation strictly reduces the number of 
tetrahedra. 

To obtain the 0(n) running time, we perform the crushing process as follows. 
We first identify which tetrahedra of T survive the process (i.e., do not contain 
quadrilaterals of S). We then identify how their faces are glued together after 
crushing, by following chains of tetrahedra through T that do not survive (i.e., 
do contain quadrilaterals of S). Each tetrahedron appears at most twice in all of 
these chains, and so the total length of all chains — and therefore the total running 
time — is bounded by 0(n). □ 

3.2. Building a one- vertex triangulation. Here we give the details of steps Q] 
and [2] of Algorithm [21 in which we build a one- vertex triangulation of K. 

As noted in the introduction, step [T] is taken directly from Hass et al. [21] , who 
show how to convert a c-crossing diagram of a knot K into a triangulation of the 
knot complement of K built from 0(c) tetrahedra in 0(c log c) time. 

What remains is to make this into a one-vertex triangulation T . There are 
well-known algorithms for producing one- vertex triangulations of a 3-manifold |27[ 
|3"3"] . though they have not been studied from a complexity viewpoint; we base our 
algorithm on subcomplex expansion and the Jaco-Rubinstein crushing technique, 
and show that its running time is a small polynomial in c. We briefly sketch the 
key ideas here, and give the full details in the appendix. 

Lemma 6. Given an n-tetrahedron triangulation T of a knot complement K , we 
can construct a one-vertex triangulation T' of K with < n tetrahedra in 0(n 3 ) 
time. 
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Figure 9. Crushing an edge between two vertices 



Proof. The key idea is to repeatedly modify the triangulation by crushing edges that 
run between two distinct vertices (see Figure [3]), until only one vertex remains. 

In a generalised triangulation it is not always possible to crush an edge e directly 
|10) (for instance, if e appears more than once in the same tetrahedron), and so 
instead we achieve this by crushing a normal surface that surrounds e using the 
Jaco-Rubinstein technique. If e is not "directly surrounded" by a normal surface 
(i.e., the frontier of a small regular neighbourhood of e is not normal), then we 
expand e to a subcomplex £ C T that is. Either way, we obtain a normal sphere 
or disc S that "cuts out" a 3-ball in T containing e. After crushing, this 3-ball 
(including e) will disappear. 

We show that building the subcomplex £ and the surface S takes 0(n 2 ) time, 
and by Lemma [5] the subsequent crushing operation is 0(n). Since each iteration 
reduces the number of tetrahedra it follows that the entire one-vertex conversion 
runs in time 0(n 3 ). See the appendix for further details. □ 

Noting that the initial triangulation of the knot complement K has 0(c) tetra- 
hedra, where c is the number of crossings in the knot diagram, we therefore have: 

Corollary 7. Steps[l]and\^ of Algorithm^ can be performed in 0(c 3 ) time, where 
the input knot diagram has c crossings, and the resulting triangulation T of K has 
0{c) tetrahedra. 

Remark. In experimentation, we did not need the complex operations outlined 
above. Instead we simplified each triangulation using the ready-made (and poly- 
nomial time) simplification tools from the software package Regina [10] . which in 
every case produced a one- vertex triangulation of K with no further work required. 

Observation 8. If T is a one-vertex triangulation of a knot complement, then 
the boundary torus dT contains precisely two triangular faces as illustrated in Fig- 
ure 



2(a) and the unique vertex of T also lies on this boundary dT. 



The proof of this well-known fact involves a simple Euler characteristic argument, 
which we do not reiterate here. 

3.3. Searching for positive Euler characteristic. We now describe step [3] of 
Algorithm [21 which sits at the heart of the algorithm (and is its major bottleneck). 

Assumptions. Throughout this section, T is a one- vertex, n-tetrahedron triangu- 
lation of a knot complement K. 

The aim of step [3] is to locate a connected normal surface S in T which is not 
the vertex link and which has positive Euler characteristic, or to show that no such 
surface exists. 

In principle, we do this by treating the Euler characteristic x : -> R as a 
linear objective function, and maximising it over the space of admissible points in 
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M. 7n using branch-and-bound techniques. If the desired surface S exists then the 
maximum \ is unbounded, and if not then the maximum x is zerod In practice, 
since we only need to distinguish between a zero or positive maximum, we formulate 
the search purely in terms of branching and feasibility tests. 

To ensure that the resulting surface S is connected, we set as many coordinates 
to zero as possible (see Lemma |9|). To ensure that S is not a vertex link, we ensure 
that at least one of its triangle coordinates is zero (see Lemma [TO]) . 

Lemma 9. Let S be a normal surface in T with positive Euler characteristic, 
and suppose there is no normal surface S' in T with positive Euler characteristic 
with the following properties: (i) if the ith coordinate of v(S) is zero then the ith 
coordinate ofv(S') is likewise zero; (ii) there is some i for which the ith coordinate 
ofv(S) is non-zero but the ith coordinate o/v(S") is zero. 

Then the smallest positive rational multiple of v(S) whose coordinates are all 
integers represents a connected normal surface with positive Euler characteristic. 

Proof. In essence, we show that such a surface S must be a vertex normal surface, 
i.e., one whose vector v(S) lies on an extreme ray of the cone in R 7 ™ defined by 
x > and the matching equations Ax = 0. We then use the fact that the smallest 
integer multiple of a vertex normal surface is always connected. See the appendix 
for details. □ 

Lemma 10. A connected normal surface S in T is the vertex link if and only if 
the corresponding vector v(S) G M 7 " has none of its An triangle coordinates equal 
to zero. 

Proof. This is an immediate consequence of the fact that T is one-vertex. The 
vertex link itself has all triangle coordinates equal to one, and any surface whose 
triangle coordinates are all positive must be the disconnected union of the vertex 
link with zero or more other surfaces. □ 

The overall structure of our search is as follows: 

Algorithm 11. To find a connected normal surface S which is not the vertex link 
and which has positive Euler characteristic: 

(1) Search for an admissible point p G M. for which x(p) > lj an d for which 
at least one of the 4n triangle coordinates is zero. If no such p exists then 
the surface S does not exist either. 

(2) Construct a system C of linear constraints, initially defined by x > 0, 
x( x ) ^ 1; an d the matching equations Ax = 0. Using the point p located 
in the previous step: 

(a) For each coordinate i — l,...,7n with pi — 0, add the additional 
constraint Xi — to the system C. 

(b) Then, for each coordinate i = 1, . . . , 7n with pi ^ 0, add the constraint 
Xi = to C and test whether C is feasible (i.e., has any solutions) . If 
so, keep the constraint Xi — in C; if not, remove it again. 

(3) Let q be a solution to the final system C, and let Aq be the smallest positive 
rational multiple of q whose coordinates are all integers. Then Aq is an 
integer vector in R 7 ™ that represents the desired surface S. 



This is because positive solutions can be arbitrarily scaled, and the zero vector is always 
admissible. 
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We give further details for step [T] shortly (this is where the branch- and-bound 
techniques are used). In the meantime: 

Lemma 12. Algorithmic is correct in either finding the desired surface S or show- 
ing that it does not exist. Moreover, steps [J2 and [3] have running time polynomial 
in n. 

Proof. Here we outline the key ideas in the proof; see the appendix for full details. 

We express positive Euler characteristic using the condition \ > 1 (not x > 0), 
because the linear programming techniques in our feasibility tests require non-strict 
inequalities. This does not affect the existence of solutions, because any admissible 
x £ M. 7n with x( x ) > can be scaled to an admissible x' = Ax 6 R 7 " with A G N 
and x(Ax') > 1. 

Step Q] uses Lemma [10] to exclude the vertex link from our search. If we do find 
an admissible point p S R 7n in step[TJ then steps [2] and [3] enforce the conditions of 
Lemma [9] to ensure that the final surface is connected. Since p is admissible, the 
additional constraints Xi = that we add in step [U[a) ensure that the final vector 
q is admissible also. 

In step[3l we prove that q does have a rational multiple with integer coordinates 
by showing that the final system C has a one-dimensional solution set. 

Steps [2] and [3] run in polynomial time because the only significant work is in 
building and solving In + 1 linear programs (the In feasibility tests in step [21 plus 
the search for a final solution in step [3]). The techniques of Hass et al. show that 
the final integer vector Aq has coordinates bounded by exp(0(n)), and so can be 
stored using 0(n) space. □ 

What remains now is to describe step Q] of Algorithm 1 1 1 1 where we search for 
an admissible point p € R 7 ™ with x(p) > 1 an d at least one of the An triangle 
coordinates equal to zero. Recall that an admissible point x g R 7n must satisfy 
x > 0, the matching equations Ax = 0, and the quadrilateral constraints. For 
convenience, let xf" denote the ith triangle coordinate of x (1 < i < 4n), and let 
x^ k denote the quadrilateral coordinate that counts the fcth type of quadrilateral 
in the ith tetrahedron (1 < t < n, 1 < k < 3). 

The conditions x > 0, As. = and x(x) > 1 can be described using a system of 
linear constraints over R . The remaining conditions require branching: 

• Triangle branching: 

We must make an initial decision on which triangle coordinate will be zero. 
To avoid redundancy, we select the first such coordinate; that is, we choose 
some i £ {1, . . . , 4n} for which x-^ , . . . , xf^-^ > 1 and xf~ = 0. 

• Quadrilateral branching: 

For each i = 1, . . . , n, we must decide which quadrilateral coordinate in the 
ith tetrahedron will be non-zero, if any. That is, we choose from (i) = 
x?,2 = = 0; (ii) x° > 1 and z° = z° = 0; (iii) z° > 1 and 
= = 0; or (iv) x% > 1 and x° = i° = 0. 

As before, we replace all strict inequalities Xi > with non-strict inequalities 
Xi > 1; again this makes the linear programming simpler, and does not affect the 
existence of solutions because any admissible solution can be rescaled to one with 
all coordinates > 1. 
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Figure 10. The combinatorial search tree 



We arrange these decisions into a combinatorial search tree, as illustrated in 
Figure [TUJ the triangle branch is chosen first, and then the n quadrilateral branches 
are chosen in some order (which we describe shortly). To reduce the total amount 
of branching, we merge quadrilateral branches (i) and (ii) into the single branch 
(i*) > and x^ 2 = x ^3 = 0, so there are only three (not four) branches at each 
quadrilateral decision. 

To run step [1] of Algorithm [TTJ we begin at the root of this search tree and 
traverse it in a depth-first fashion. At all stages we maintain a system Af of linear 
constraints corresponding to the branches that have been chosen: 

• At the root of the tree we initialise Af to the constraints x > 0, the matching 
equations Ax. = 0, and x( x ) > 1- 

• Each time we follow a new branch down, we extend Af by adding addi- 
tional constraints to describe the branch that was chosen. For instance, 
when following a triangle branch we would add constraints of the form 
x 1 , . . . , > 1 and xf" = 0, and when following the third quadrilateral 
branch for some tetrahedron we would add the constraints x^ 3 > 1 and 

• Conversely, each time we backtrack and follow a branch up, we remove 
these additional constraints from Af. 

Crucially, each time we enter a node in the search tree we test the system Af 
for feasibility (i.e., whether it has a solution). If the system is not feasible, then 
we backtrack immediately (since the branches chosen thus far are inconsistent, and 
cannot lead to the desired admissible point p £ M. 7n ). If we ever reach a leaf 
node of the tree — that is, a node at which we have chosen a triangle branch plus 
quadrilateral branches for all n tetrahedra — then any solution to the system Af is 
the admissible point p G R 7 ™ that we require. 

It is clear from the layout of the search tree and the structure of the system 
Af that this method is correct in performing step [T] of Algorithm 1111 Regarding 
running time: 

Lemma 13. Using the branching method described above, step [71 of Alaorithm \ll\ 
runs in time 0(3™ x poly(n)). 

Proof. We simply count nodes: there are 4n branches for the initial triangle decision 
and three branches for each of the n quadrilateral decisions, giving 4n-3" leaf nodes 
in the worst case. Each feasibility test can be run in polynomial time using linear 
programming techniques, yielding an overall running time of 0(3™ x poly(n)). □ 
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Combining this result with Lemma IT2| we have: 

Corollary 14. Algorithm ] 11\ [i.e., step\3\ of Algorithm^! is correct in either finding 
the desired surface S or showing that it does not exist, and runs in time 0(3™ x 
poly(n)). 

Again, we highlight the underlying branch-and-bound motivations. If we treat 
X as a linear objective function to maximise, then our feasibility tests at each node 
play the role of linear relaxations — in the absence of constraints for decisions not yet 
made, they test whether the maximum \ could possibly be positive in the subtree 
below this node. 

We finish with important enhancements to the vertical layout of the search tree, 
i.e., the order in which we make our n quadrilateral decisions. Figure [10] shows these 
decisions ordered by tetrahedron number. In reality, we dynamically reorder these 
decisions to take full advantage of the interaction between different constraints and 
to prune the search tree as frequently as possible. Experimentation suggests that 
these optimisations are crucial to the polynomial-time behaviour that we observe 
in Section [5] 

• After choosing an initial triangle branch, if the triangle coordinate that 
we set to zero was in the zth tetrahedron, then we immediately branch on 
quadrilateral coordinates in the same (ith) tetrahedron. 

This is because setting a triangle and two quadrilateral coordinates to 
zero in the same tetrahedron is an extremely powerful constraint, and im- 
mediately eliminates several other "nearby" quadrilateral and triangle co- 
ordinates from surrounding tetrahedra. 

• For later quadrilateral decisions, we greedily branch on the tetrahedron 
that yields the fewest feasible child nodes. That is, for each unused tetra- 
hedron we build the constraint systems Af that would result at each of the 
three child nodes if we branched next on that tetrahedron, and count how 
many of these three systems are feasible. We then branch on quadrilateral 
coordinates in the tetrahedron that minimises this count. 

Essentially, this greedy approach allows us to lock in additional "forced" 
constraints as quickly as possible. Moreover, if there is some tetrahedron 
in which none of the three quadrilateral branches yield a feasible child 
node, our approach detects this and allows us to backtrack immediately. 
A drawback of this greedy method is that it requires a linear number of 
feasibility tests at each node of the search tree, but this only affects the 
running time by a polynomial factor. 

3.4. Crushing normal surfaces. Here we finalise the details for step [4] of Algo- 
rithm d 

Corollary 15. All claims made in step^ of Algorithm^ are correct, and this step 
can be performed in 0(n) time. 

Proof. Regarding correctness: By Theorem [TJ it is clear that if the surface S does 
not exist then K is non-trivial, and if S is a disc with non-trivial boundary then K 
is trivial. Otherwise, we crush S: here Lemma [5] shows that the new triangulation 
7~~' represents the same knot complement K with strictly fewer than n tetrahedra. 

Regarding running time: Since S (if it exists) is connected with positive Euler 
characteristic, and since T is a one- vertex triangulation whose torus boundary has 
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only two triangles (see Observation [8]), then S is a disc with non-trivial boundary 
if and only if S meets dT in any non-empty curve other than the trivial curve 
that encircles the vertex. This is easily tested in 0(n) time just by examining the 
coordinates of v(5). Lemma [5] shows that the crushing can likewise be performed 
in 0(n) time. □ 

3.5. Proofs of correctness and running time. To prove that the full unknot 
recognition algorithm is correct and has the claimed running time, all that remains 
now is to tie together results from the previous sections. 

Proof of Theorem^ (correctness): This follows immediately from Corollaries 1141 
and US □ 

Proof of Theorem^ (running time): Let c be the number of crossings in the input 
knot diagram. Corollary [7] shows that steps [T] and [2] run in time polynomial in c, 
and that the resulting number of tetrahedra n is O(c). Corollary 1151 shows that 
step S] runs in polynomial time, and Corollary [14] shows that step [3] runs in time 
0(3" x poly(n)). Since each iteration of step |4] either terminates the algorithm or 
reduces the number of tetrahedra, we loop through steps [2HH at most n G 0(c) 
times. □ 



4. 3-SPHERE RECOGNITION AND PRIME DECOMPOSITION 

We can adapt our new algorithm for use with other topological problems, such 
as 3-sphere recognition and prime decomposition. We briefly sketch the key ideas 
here. 

Given a triangulation T that represents a closed orientable 3-manifold, 3-sphere 
recognition asks whether this manifold is topologically equivalent to the 3-sphere, 
and prime decomposition breaks this 3-manifold into a collection of "prime factors" 
(which are combined together using the topological operation of connected sum). 

For both problems, early algorithms were mathematical breakthroughs but al- 
gorithmically cumbersome [28, 37]. They have since enjoyed great improvements in 
both implementability and efficiency, though like unknot recognition they still run 
in worst-case exponential time. See |10j for a modern formulation of the algorithms 
as they appear today. 

For both algorithms, the central operations — and their exponential-time bottlenecks — 
are steps that search for normal or "octagonal almost normal" spheres within a tri- 
angulation T ■ An octagonal almost normal surface is like a normal surface, but in 
additional to triangles and quadrilaterals we require precisely one octagonal piece 
in one tetrahedron. 

We can locate normal spheres directly using Algorithm ITTI four new search based 
on branching and feasibility tests), and to locate almost normal spheres we use a 
variant that is modified to work with almost normal coordinate systems. These 
new search algorithms can be dropped directly into the modern 3-sphere and prime 
decomposition algorithms, replacing traditional search techniques based on vertex 
enumeration in high-dimensional polytopes. 

5. Experimental performance 

Here we describe the results of extensive testing of the new unknot recognition 
algorithm over all 2977 prime knots with < 12 crossings. We also test the new 
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algorithm for 3-sphere recognition, as outlined in Section [4] by running it over the 
first 150 homology spheres in the Hodgson- Weeks census [24j . 

The algorithms have been implemented in C++ in the open-source computational 
topology software package Regina [71 111], and will be included as part of Regina's 
next public release (version 4.94, due in early 2013). We briefly describe some 
aspects of the implementation, and then describe the experimental data and results. 

5.1. Implementation. Although it is possible to test the feasibility of a system 
of linear constraints in polynomial time (for instance, using interior point methods 
[29] , we use a variant of the simplex method due to its ease of implementation and 
its excellent performance in practical settings (as discussed further in Section [6]) . 
Specifically, we use the revised dual simplex method |31j with an implementation 
that exploits the sparseness of the matching equations Ax = 0. 

For a pivoting rule, we use Dantzig's classical method of choosing the exiting 
variable with largest magnitude negative value in the tableaux [15]. Although this 
method is fast and works well to reduce the total number of pivots, it can lead to 
cycling. We therefore use Brent's algorithm to detect cycling [5], and when it occurs 
we switch to Bland's rule instead [4], which exhibits weaker performance but does 
not cycle. We note that cycling was indeed detected during our experimentation, 
and so these cycle-breaking techniques are necessary in practice. 

All computations use exact integer and rational arithmetic, provided by the GNU 
multiple precision arithmetic library [19] . To limit the overhead, we work in native 
integers wherever possible but test for overflow on all arithmetical operations, and 
only switch to exact arithmetic when necessary. This behaviour, which improves 
performance surprisingly well, is inspired by (but far less sophisticated than) the 
lazy evaluation methods used for exact arithmetic in the CGAL computational 
geometry library [6j Q] . 

5.2. Experimental results. For unknot recognition, our experimental data set 
consists of all 2977 prime knots with < 12 crossings, as taken from the Knotlnfo 
database. For 3-sphere recognition, our data set consists of all 150 homology spheres 
from the Hodgson- Weeks census of closed hyperbolic 3-manifolds [24] . Both are 
intended as exhaustive and "punishing" data sets, where simplification tools cannot 
solve the problems (since the inputs are not trivial knots or 3-spheres), and where 
our branching search needs to conclusively determine that certain normal or almost 
normal surfaces (such as a disc with non-trivial boundary) do not exist. 

As is common in computational topology, our first step before running any other 
algorithms is to simplify the input triangulations using a suite of local moves; we 
use the ready-made (and polynomial time) suite from Regina [10] . The resulting 
triangulations are large: the knot complements have up to n = 50 tetrahedra, and 
the homology spheres have up to n = 20 tetrahedra. It is a testament to the 
strength of the simplification suite that for all 2977 + 150 inputs we only require a 
single pass through Algorithm [TTJ (the branching search) — we never need to crush 
away "junk" discs or spheres and run the search again (as in step [4] of Algorithm [2]) . 

Figure [TT] counts how many nodes the new algorithms visit in the search tree (as 
described in Section 13.31) H this measure is crucial because it determines how many 



What "visit" means in this context relies on details of the implementation, but up to asymp- 
totics this is unimportant. 
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Figure 1 1 . The number of nodes visited in the search tree 
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Figure 12. Summary of running times for unknot recognition and 
3-sphere recognition 

feasibility tests we run, and is the source of the theoretical exponential running 
time. For 3-sphere recognition the results are noisy, but for unknot recognition 
they are unequivocally linear: in every case the number of nodes lies between 8n 
and lOrt, and for all but 24 of the inputs the figure is precisely 8n + 1 (essentially 
indicating that we never branch on quadrilateral coordinates at all). This linear 
growth in the number of nodes corresponds to a polynomial running time, and 
explains the exceptional performance of the algorithm. 

We measure the running times of the new algorithms against prior state-of- 
the-art algorithms based on normal surface enumeration, including (i) traditional 
algorithms based on the double description method [5], and (ii) recent algorithms 
based on tree traversal |12) . These prior algorithms are also implemented in Regina, 
and their code is heavily optimised; see |10) for an overview of each. 

Figurell"2"lsummarises the running times (single-threaded, measured on a 2.93 GHz 
Intel Core i7); note that the time axis is logarithmic. For each number of tetrahedra 
n we aggregate times using geometric mean (which benefits the prior algorithms, 
which have much wider variability) . See Figures [T3] and [TJ] in the appendix for 
detailed plots with individual data points. 

The results are extremely pleasing. For 3-sphere recognition, the new algorithm 
runs over 15 x faster than the prior tree traversal algorithm and roughly 223 x faster 
than the prior double description algorithm for the final case (n = 20). Moreover, 
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as seen in the appendix, the new algorithm exhibits much less variability in running 
times. 

However, it is with unknot recognition that the new algorithm truly shines. Here 
we only ran the prior algorithms over the first 515 input knots (as sorted by increas- 
ing n). This is because running times became too slow to proceed any further — a 
linear regression on log (time) suggests that for the largest case with n = 50 the 
prior tree traversal algorithm would take ~ 6000 years, and the prior double de- 
scription method would take ~ 40, 000 years. In contrast, the new algorithm solved 
all 2977 input cases (including the largest case with n — 50) in under 5 minutes 
each. 

6. Discussion 

Although our new algorithm for unknot recognition remains exponential time 
in theory, the observed polynomial growth in practice is extremely pleasing, and 
indeed very exciting — this moves the study of unknot recognition into a new phase, 
where we can now solve it conclusively and quickly "in practice" , albeit without the- 
oretical guarantees on the running time. This is reminiscent of the simplex method 
for linear programming, a fast-in-practice but exponential-in-theory algorithm for 
linear programming which was known long before polynomial alternatives were dis- 
covered [201 [29] , and which still enjoys widespread use because of its extremely fast 
"typical" behaviour. 

We have a good understanding of why the simplex method works well in practice: 
it has since been shown to be polynomial time in settings such as average, generic 
and smoothed complexity (38j [39] [40] . In contrast, there are no such results for 
unknot recognition; more generally, results of this type are extremely scarce in 
the study of topological algorithms on 3-dimensional triangulations (the setting for 
many other knot algorithms also). Reasons for this include: 

• Random 3-manifold triangulations are difficult to study because they are 
extremely rare: the probability that a random pairwise gluing of faces from 
n tetrahedra yields a 3-manifold triangulation approaches zero as n — > oo 
[17] . and there is no known algorithm for generating a random 3-manifold 
triangulation (under some useful probability distribution) in polynomial 
time [36] . 

• "Walking" through the space of 3-manifold triangulations is difficult be- 
cause the diameter of this space could be extremely large: for knot com- 
plements the best known bounds involve exponentially high towers of ex- 
ponentials [54] , 

In Section [4] we adapted the new unknot recognition algorithm to the problems 
of 3-sphere recognition and prime decomposition. These techniques have potential 
to extend to an even broader range of topological problems: promising candidates 
include 0-efficiency algorithms, and the difficult but important problem of finding 
incompressible surfaces. 
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Appendix A. Detailed experimental data 

In this first appendix we present the individual case-by-case experimental data 
points that were aggregated to form the plots in Section [SJ 

Figure [T3] presents the individual running times for unknot recognition over all 
2977 prime knots with < 12 crossings. The upper two plots show these times for the 
two prior state-of-the-art algorithms (one based on the double description method 
[8j, and one based on tree traversal [32]). The lower left plot shows running times 
for the new Algorithm [2l Every point in these plots represents a single input knot 
K. 

As explained in Section [5l the prior algorithms (the upper two plots) only show 
times for the first 515 input knots, since running times quickly reached the point 
of infeasibility. The new algorithm shows times for all 2977 inputs. 

The lower right plot aggregates times for each n by geometric mean: this is 
identical to the summary in Figure [T^] from the main body of this paper, and is 
included again here for convenience. 

In a similar fashion, Figure [14] presents individual running times for 3-sphere 
recognition over all 150 input homology spheres (this time all three algorithms 
show results for all 150 inputs). 

Appendix B. Additional proofs 

Here we give full details for several proofs that were only sketched briefly in the 
main body of the paper. These proofs are for Lemmas l5l \6\ l9l and [T2l 

Lemma 5. Given a triangulation T of a knot complement K with n > 1 tetrahedra 
and a connected normal surface S in T which is not a vertex link and which has 
positive Euler characteristic, we can crush S using the Jaco- Rubinstein technique 
and convert the result into a new triangulation T' of K with strictly fewer than n 
tetrahedra. Moreover, this process can be carried out in 0(n) time. 
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Figure 13. Running times for unknot recognition on the first 2977 
prime knots 



Proof. The only connected surfaces with positive Euler characteristic are the sphere, 
the disc, and the projective plane. Since T triangulates a knot complement in S 3 
and S 3 does not contain an embedded projective plane, we conclude that S is either 
a sphere or a disc. 

By results of Jaco and Rubinstein |27j . crushing a normal sphere or a disc in 
an orientable triangulation with boundary results in a new triangulation (possibly 
disconnected or possibly empty) whose underlying manifold is obtained from K by 
zero or more of the following operations: 

(i) cutting manifolds open along spheres and filling the resulting boundary spheres 
with 3-balls; 

(ii) cutting manifolds open along properly embedded discs; 
(hi) capping boundary spheres of manifolds with 3-balls; 

(iv) deleting entire connected components that are any of the 3-ball, the 3-sphere, 
projective space MP 3 , the lens space L 3i i or the product space S 2 x S 1 . 

Since K is a knot complement, we can make this list significantly simpler: 

(i) Any embedded sphere in a knot complement bounds a ball, and so the first 
operation is equivalent to adding a new (disconnected) 3-sphere component. 

(ii) If D is a properly embedded disc whose boundary is trivial on dK, then 
since K is a knot complement it must be true that D and the corresponding 
section of dK together bound a 3-ball, and the operation of cutting along D 
is equivalent to adding a new 3-ball component. 
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Figure 14. Running times for 3-sphere recognition on the first 
150 homology spheres 



If D is a properly embedded disc whose boundary is non-trivial on dK 
then K must be a trivial knot, and cutting along D converts the manifold K 
into a 3-ball. 

(iii) The only boundary spheres that might occur will be the boundaries of new 
3-ball components created in the previous step. Capping these spheres with 3- 
balls has the effect of converting our new 3-ball components into new 3-sphere 
components instead. 

(iv) If we begin with a knot complement, none of the operations above can ever 
produce the manifolds RP 3 , L3.1 or S 2 x 5 1 . Therefore the only components 
that we might ever delete are 3-balls and 3-spheres. 

In other words, the topological effects of crushing T along the sphere or disc S 
are limited to replacing K with a 3-ball (but only if the knot K is trivial), and 
adding or removing 3-ball and/or 3-sphere components. 

It is clear now how to obtain the new triangulation T 1 of K. If, after crushing, 
any connected component of the resulting triangulation has torus boundary then we 
can take this as the new triangulation T' ■ Otherwise we know that K is trivial, and 
we take T' to be the standard one-tetrahedron triangulation of the solid torus [27] 
(though of course we could just as well terminate the unknot recognition algorithm 
immediately) . Either way, we have that T 7 triangulates the same knot complement 
K. 
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Because S is not the vertex link, it contains at least one quadrilateral and so 
crushing will strictly reduce the number of tetrahedra. It is clear then that after 
extracting the component with torus boundary (or building a new one-tetrahedron 
solid torus), the number of tetrahedra in T' will be strictly less than n. 

It remains to show that this procedure can be carried out in 0(n) time. There 
are two stages: (i) the Jaco-Rubinstein crushing process, and (ii) extracting the 
connected component with torus boundary (if any). We focus first on stage (ii), 
where the 0(n) time is simple to obtain through a series of standard amortised 
depth-first search arguments: 

• Since each tetrahedron is adjacent to at most four neighbours, we can iden- 
tify connected components in 0(n) time by following a depth- first search 
through adjacent tetrahedra for each component. This takes 0(t) time for 
each i-tetrahedron component, summing to 0(n) time overall. 

• For each connected component, we "build the skeleton"; that is, group the 
At vertices, 6t edges and At faces of the t individual tetrahedra into equiv- 
alence classes that indicate how these are identified (or "glued together") 
in the triangulation. As before, we do this by following a depth-first search 
through adjacent tetrahedra for each equivalence class. Each search re- 
quires time proportional to the size of the equivalence class, summing to 
0(t) time for each component and 0(n) time overall. 

• For components with non-empty boundary, we can test for torus boundary 
by counting equivalence classes of vertices, edges and faces, and computing 
the Euler characteristic \ — vertices — edges + faces — tetrahedra. This is 
enough to distinguish between components with torus boundary (x = 0) 
and sphere boundary (x = l)j which from earlier are the only possible 
scenarios after crushing. Once again this sums to 0(n) time overall. 

Our remaining task is to show that the Jaco-Rubinstein crushing process can be 
carried out in 0(n) time. It is simple to identify which tetrahedra of T survive the 
crushing process (those whose corresponding quadrilateral coordinates in S are all 
zero) — the challenge is to identify how the faces of these surviving tetrahedra are 
to be glued together. 

Let <fi be some face of a surviving tetrahedron A. To identify the new partner 
face for 0, we follow a path through adjacent tetrahedra in the original triangula- 
tion T, as illustrated in Figure [151 Whenever we enter a tetrahedron that contains 
quadrilaterals of S, we cross to the face on the opposite side of these quadrilaterals 
and continue through to the next adjacent tetrahedron. If we ever reach a tetrahe- 
dron with no quadrilaterals of S, then the resulting face <j>' is the partner to which 
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4> is glued. If instead we reach the boundary 97", then becomes a boundary face 
of the new triangulation. 

To see why following paths in this way correctly identifies the face gluings, we 
must study the way in which pillows (formed on either side of the quadrilaterals in 
a tetrahedron) flatten to triangles, as seen earlier in Figure [8] Jaco and Rubinstein 
explain this in their original paper [2 7) . and we do not reiterate the details here. It 
is important to note that such a path cannot cycle (since the starting tetrahedron 
A does not contain quadrilaterals) , and so we must eventually reach one of the two 
conclusions described above. 

Returning to time complexity: although any individual path may be long, we 
note that each tetrahedron of T can feature at most twice amongst all the paths 
(once for each of the two directions in which we might cross over the quadrilaterals) . 
Therefore the total length of all such paths is 0{n). Each step in such a path takes 
0(1) time to compute (since all we need to know is which, if any, of the three 
quadrilateral coordinates is non-zero in the next tetrahedron), and so the gluings 
in the final crushed triangulation can be computed in total O(n) time. □ 

Lemma 6. Given an n-tetrahedron triangulation T of a knot complement K , we 
can construct a one-vertex triangulation T' of K with < n tetrahedra in 0(n 3 ) 
time. 

Proof. First, note that in 0(n) time we can partition the 4n individual tetrahedron 
vertices into equivalence classes that indicate how they are identified in the overall 
triangulation, using the same depth-first search through adjacent tetrahedra as seen 
earlier in the proof of Lemma [5] As a result, in 0(n) time we can test whether T 
has more than one vertex and, if it does, locate an edge e that joins two distinct 
vertices of T. 

Suppose then that T has more than one vertex, and let e be an edge that joins 
two distinct vertices in T. If T contains more than one vertex on the boundary 
then we choose e to lie entirely in dT (as opposed to cutting through the interior 
of T); otherwise we choose an edge e that joins the boundary vertex with some 
internal vertex (as opposed to joining two internal vertices). 

We observe that T must have n > 2 tetrahedra. This follows from a simple 
enumeration of all one-tetrahedron triangulations: the only one-tetrahedron knot 
complement is the standard one-tetrahedron triangulation of the solid torus [27] . 
which has just one vertex. 

Our plan is to convert this edge e into a normal disc S that is not a vertex link, 
all within 0(n 2 ) time. Given such a surface, Lemma [5] shows that we can crush it 
in 0(n) time to obtain a new triangulation of K with strictly fewer tetrahedra. We 
repeat this process until we arrive at a one-vertex triangulation; since the number 
of tetrahedra decreases at every stage, the process must terminate after at most n 
iterations, giving a running time of 0(n 3 ) overall. 

The remainder of this proof is devoted to showing how to build this normal disc 
S. Ideally we would take S to be the surface "immediately surrounding" e (that 
is, the frontier of a regular neighbourhood of e); then crushing S would have the 
effect of crushing the edge e itself, as illustrated earlier in Figure [9] Unfortunately 
such a surface S need not be normal — in particular, if e appears as two edges of 
some triangular face of T, then S will meet the face in a "bent" arc that does not 
correspond to any normal disc, as shown in Figure [1^1 
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Figure 16. A non- normal surface surrounding e 





Figure 17. Expanding the subcomplex £ 



Our strategy then is to expand e to a subcomplex £ of T (formed from a union 
of edges, faces and/or tetrahedra of T) that is surrounded by a normal surface. To 
do this, we initialise £ to the single edge e and then repeating the following steps: 

• If any triangular face has cither two or all three of its edges in £, we expand 
£ to include the entire face (as illustrated in Figure [TTJ) ; 

• If any tetrahedron has all four of its faces in £ , we expand £ to include the 
entire tetrahedron. 

It is clear that this expansion can be carried out in 0(n 2 ) time: each individual 
expansion step can be tested and performed in 0(n) time, and since T contains 
0(n) edges, faces and tetrahedra, there can be at most 0{n) such expansion steps 
in total. 

From here on we use the phrase u £ is surrounded by S" to mean (as usual) that 
the surface S is the frontier of a small regular neighbourhood of £. If S is a disc, 
then we use the phrase "S cuts out a 3-ball" to indicate that cutting T along S 
has the effect of slicing off a 3-ball component (in other words, S bounds a 3-ball 
together with some disc in dT). 

We prove by induction that, at each stage of the construction, £ is surrounded 
by a surface S which is a disjoint union of spheres and/or discs, at least one of 
which is a disc that cuts out a 3-ball containing e. 

• At the beginning of the construction this is true because e joins two distinct 
vertices, at least one of which is on <9T, and because we chose e so that if 
both vertices lie on dT then so does the entire edge e. 

• If a triangular face has two edges in £ , the resulting expansion step simply 
"pushes" the surface S across the face, as shown in Figure 18(a) This does 
not introduce any new intersections with e, and therefore preserves all of 
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(a) Pushing S across a face 




(b) Cutting 5 through a face 



Figure 18. How expanding £ affects the surface S 




FIGURE 19. The effects of cutting a tube in S 

the required properties. If the face meets the boundary then we might split 
a disc into two discs, but again the necessary properties are preserved. 

• If a triangular face has three edges in £, the expansion has the effect of 
cutting a tube in S, as shown in Figure |T8(b) | This either replaces a sphere 
of S with two spheres, or replaces a disc of S with a disc and a sphere. 

If this surgery takes place on the disc that cuts out a 3-ball containing 
e, then the resulting disc will again cut out a 3-ball containing e. This is 
because (depending on whether the ball is on the inside or the outside of 
the tube) either the surgery carves off a piece of the ball (but not the piece 
that meets the boundary, and hence not the piece containing e), or else it 
attaches the contents of some embedded sphere in K (which, since K is 
a knot complement and the attached contents do not meet the boundary, 
must be another 3-ball). See Figure [19] for illustrations. 

• If a tetrahedron has four faces in £, then the resulting expansion simply 
removes a sphere from the collection. 

We now claim that, after the expansion is finished, we can take S — the frontier 
of a small regular neighbourhood of £ — to be a normal surface. Indeed, we can 
explicitly construct this surface from £ by inserting the following triangles and 
quadrilaterals into each tetrahedron A: 
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Figure 20. Building a normal surface surrounding £ 

• If £ contains the tetrahedron A, there are no normal triangles or quadri- 
laterals in A at all. 

• Otherwise, £ contains at most one face of A (since two or more faces would 
cause the expansion process to eventually consume A completely). If £ 
contains such a face, then we place a normal triangle beside this face. 

• If £ does not contain any faces of A, it might still contain either one edge 
of A, or two opposite edges of A (any more would again cause further 
expansion). Iff contains such edge(s), then we place a normal quadrilateral 
beside each edge. 

• If £ contains any additional vertices of A that are not yet accounted for, 
we place a triangle next to each such vertex. 

See Figure HOI for some illustrations of this process. It is clear that we can 
construct this surface S in 0(n) time, since it involves O(l) operations in each 
tetrahedron. Moreover, since this surface contains 0(n) triangles or quadrilater- 
als in total (at most four in each tetrahedron), we can in 0(n) time extract the 
connected components of «S, compute their Euler characteristics, and identify one 
such component that is a disc but not a vertex link. We know there is at least one 
such component to be found: the disc that cuts out a 3-ball containing e cannot 
be a vertex link, since otherwise it would cut out a 3-ball on both sides (one side 
containing e, the other containing the linked vertex), whereupon K would simply 
be a 3-ball, not a knot complement. 

This concludes our 0(n 2 ) time construction of a normal disc in T that is not a 
vertex link, and thereby concludes the proof. □ 

Lemma 9. Let S be a normal surface in T with positive Euler characteristic, 
and suppose there is no normal surface S' in T with positive Euler characteristic 
with the following properties: (i) if the ith coordinate of v(S) is zero then the ith 
coordinate ofv(S') is likewise zero; (ii) there is some i for which the ith coordinate 
ofv(S) is non-zero but the ith coordinate o/v(S") is zero. 

Then the smallest positive rational multiple of v(S*) whose coordinates are all 
integers represents a connected normal surface with positive Euler characteristic. 

Proof. Recall that a point x G M 7n is admissible if it satisfies x > 0, the matching 
equations Ax — 0, and the quadrilateral constraints. Let C denote the polyhedral 
cone in R 7n defined by x > and Ax = 0. We first show that the vector v(S') 
lies on an extreme ray of C (in the language of normal surface theory, such an S is 
called a vertex normal surface). 

Suppose that v(S) does not lie on an extreme ray of C. Since v(5) is admissible 
we still have v(S*) £ C, and so we can express v(S) as a finite non-negative linear 
combination v(5) = Aiei + . . . + Afce^, where each ei lies on a different extreme 
ray of C, each A^ > 0, and there are k > 2 terms. Because C is a rational cone, 
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we can take each ej to be an integer vector. Since C lies in the non- negative 
orthant, it follows that whenever the ith coordinate of v(S) is zero then the ith 
coordinate of each must likewise be zero. In particular, since v(S) satisfies the 
quadrilateral constraints, each also satisfies the quadrilateral constraints, and so 
each = v(S'i) for some normal surface Si. 

Recall that the Euler characteristic x is a homogeneous linear function on R 7n . 
Since x(v(S 1 )) > 0, we have x( e p) > for some summand e p ; i.e., the corresponding 
surface S p has positive Euler characteristic. Take any other summand e q where 
q 7^ p. Since the inequalities that define C are all of the form Xi > 0, the extreme 
rays of C are defined by which coordinates are zero and which are non-zero; in 
particular, there must be some coordinate which is non-zero in e q but zero in e p . 
This coordinate is therefore non-zero in v(S) but zero in v(S , p ), and so the surface S p 
satisfies all of the properties of S' in the lemma statement, yielding a contradiction. 

Therefore v(S) lies on an extreme ray of C. Let u be the smallest positive rational 
multiple of v(5) with integer coordinates. It is clear that u is likewise admissible 
with x(u) > 0; that is, u represents some normal surface U with positive Euler 
characteristic. All that remains to show is that U is connected. 

If U is not connected, we can express it as the disjoint union of non-empty normal 
surfaces U = X U Y, whereupon u = v(X) + v(Y). Since u lies on an extreme ray 
of C, it must be that both v(X) and v(Y) are smaller integer multiples of u, in 
contradiction with the definition of u. □ 

Lemma 12. Algorithm[T^is correct in either finding the desired surface S or show- 
ing that it does not exist. Moreover, steps [H and [3] have running time polynomial 
in n. 

Proof. Recall that, in step Q] of the algorithm, we search for an admissible point 
p £ R 7ra for which x(p) > L an d for which at least one of the 4n triangle coordinates 
is zero. 

If there is some connected normal surface S which is not the vertex link and which 
has positive Euler characteristic, then (using Lemma [T0|) the vector v(S) £ M. 7n is 
such a point. Therefore, if we fail to find such a point p in step[TJ the algorithm is 
correct in concluding that such a surface S does not exist. 

Suppose then that we do find such a point p in step Q] Consider now the point 
q that we construct in steps [5] and [3] of the algorithm. It is clear that this point 
exists, since the initial system C has p as a solution, each additional constraint 
added in step EJa) is again satisfied by p, and then each additional constraint 
added in step[2]Jb) is explicitly tested for feasibility. 

We claim that q is admissible. Certainly it satisfies q > and the matching 
equations A(\ — 0, since these are built into the system C. Moreover, by stepUJa) 
we know that each coordinate that is zero in p is also zero in q, and so because 
p satisfies the quadrilateral constraints then q must also. Therefore q is indeed 
admissible. 

We show now that q scales to a rational vector. Let V be the rational polyhedron 
in R 7n defined by the system L. From the structure of the inequalities that define 
C, we see that each facet of V is obtained by intersecting V with a supporting 
hyperplane of the form (i) xi = 0, or (ii) x( x ) — 1- From step [2] we see that any 
intersection of V with a hyperplane Xi = either includes all of V or is the empty 
set; either way we cannot obtain a facet. Therefore V has only one facet (the 
intersection with x( x ) = 1); an d it follows that V is a one-dimensional ray. Since 
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£ is a rational system we conclude that the solution q scales down to a rational 
point (the vertex of V at the beginning of this ray). 

It follows that in step [3] the multiple Aq is well-defined, and represents an admis- 
sible integer vector with x(Aq) > 0- We therefore have Aq = v(S) for some normal 
surface S in T with positive Euler characteristic. Because some triangle coordinate 
of p is zero, stepEJa) ensures that some triangle coordinate of v(5) is zero, and it 
follows from Lemma [TOl that S is not the vertex link. 

We now use Lemma to show that this surface is connected. Suppose there 
is some normal surface S' with positive Euler characteristic where (i) for every 
coordinate of Aq = v(5) which is zero, the corresponding coordinate of v(S") is 
zero; and (ii) there is some i for which the ith coordinate of Aq = v(S) is non-zero 
but the ith coordinate of v(5") is zero. Then in step [2jb) of the algorithm, for- 
tius value of i, the constraint Xi — would have given a feasible system (with 
v(5") as a possible solution). Therefore Xj = would appear in the final system C, 
contradicting the fact that q is a solution to L. Therefore S satisfies the conditions 
of Lemma El and so (by construction of the scaling factor A) S is connected. 

We now have that S is a connected normal surface which is not the vertex link 
and which has positive Euler characteristic, and so the output of the algorithm is 
correct. All that remains is to show that steps [5] and [3] run in polynomial time. 

Step [5] runs in polynomial time because throughout its evolution the system 
C always contains 0{n) constraints, we only modify it O(n) times, and for each 
modification we run a single feasibility test which can be performed using standard 
polynomial time linear programming techniques. 

Because the final system C is one-dimensional, and because x( x ) > 1 is the only 
non-homogeneous constraint in £, it follows that removing the constraint x( x ) > 1 
has the effect of extending the solution set to a ray in R 7n that emanates from the 
origin. If we then remove the constraints n — that were inserted in step [21 the 
solution set enlarges to a polyhedral cone (defined by x > and Ax = 0) for which 
our solution q lies on an extreme ray. 

Hass et al. |23l Lemma 6.1] show that any extreme ray of this cone can be 
expressed as an integer vector with 0{n) bits — that is, with coordinates bounded 
by exp(0(n)) — and so we can solve for q in step [3J and perform the subsequent 
arithmetic in polynomial time. □ 

School of Mathematics and Physics, The University of Queensland, Brisbane QLD 
4072, Australia 

E-mail address: babamaths.uq.edu.au 

School of Mathematical and Geospatial Sciences, RMIT University, GPO Box 2476V, 
Melbourne VIC 3001, Australia 

E-mail address: melih.ozlenarmit.edu.au 



